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ABSTRACT 



Context. The radial structure of accretion disks is a fundamental issue regarding star and planet formation. Many theoretical studies, 
focussing on different aspects such as e.g. disk emissivity or ionization, have been conducted in the context of the Standard Accretion 
Disk (SAD) model, where no jet is present. 

Aims. We wish to calculate the structure of YSO accretion disks in an approach that takes into account the presence of the protostellar 

jets. The radial structure of these Jet Emitting Disks (JED) should then be compared to that of standard accretion disks. 

Methods. The analytical treatment used in this work is very similar to that of standard accretion disks but is using the parameter 

space of Magnetised Accretion-Ejection Structures that include the jet torque on the underlying disk. In this framework, the analytical 

expressions of key quantities, such as mid-plane temperatures, surface densities or disk aspect ratio are derived. 

Results. It is found that JEDs present a structure very different from the S ADs and that can be observationally tested. The implications 

on planet formation in the inner regions of accretion disks are briefly discussed. We also supply sets of analytical formulae, valid in 

different opacity regimes, for the disk quantities. These expressions can be readily used for any work where the disk structure is 

needed as an input for the model. 

Key words, accretion, accretion disks — ISM: jets and outflows — Stars: formation - method: analytical 



1. Introduction 

Accretion disks are ubiquitous in the Universe. In particular, 
they are found in Active Galactic Nuclei (AGN), around stel- 
lar black holes (X-ray binaries) and in Young Stellar Objects 
(YSO). From the theoretical point of view, accretion disks 
have been extensively studied in the c ontext of the Standard 
Accretion Disk model (hereafter SAD , iPringle & ReesI 119721: 
IShakura & SunvaevI 19731: iPrindd 19811) . be it for AGN or YSO. 
These early studies were conducted before the discovery of jets 
and focused on the physics of accretion only. Many refinements 
have been included to the initial SAD approach over the last 
thirty years, but the basic idea is still that of the first seminal 
papers: the gas inward motion is ensured by the radial turbulent 
transport of angular momentum from the inner to the outer parts 
of the disk. T o date, the magneto-rotational instability (MRI, 
iBalbus & Haw lev 1991) is the best candidate to provide and sus- 
tain the level of turbulence required in accretion disks. For the 
specific case of YSO — which is our main concern in this work — 
the standard theory of accretion disk has been widely used to de- 
termine the radial and vertical structure of accretion disks. This 
is a fundamental issue if one is to understand how star form but 
also the initial conditions of planet formation and migration. 

However, after a few decades of observations, it has ap- 
peared very clearly that disk accretion onto a central object and 
bipolar ejections cannot be disentangled. Very briefly, accretion 
is believed to power the jets which, in turn, vertically remove 
part of the disk angular momentum allowing accretion to pro- 
ceed. This accretion-ejection picture is observed on many as- 
trophysical scales: jets are present in AGN, microquasars, YSO 
and hav e more recen t ly been observed emerging f rom brown 
dwarfs (lFerrarilll998t Mirabel & Rodrigue j [l999t iBallv et al.l 
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120071: IWhelan et al.ll2005h . Despite the advances in the standard 
theory of accretion disk, the latter does not provide any explana- 
tion to the production of jets. 

In this work, we will focus on the specific case of YSO. 
About 30% of T Tauri stars (Class 2 objects) present bipo- 
lar ejection. This percentage increases to 100% for Class 
objects, the earliest stage of star formation. Several models 
have been developed to explain th e jets seen in T T auri stars. 
Stellar winds have been invoked (ISauty et al.ll2002h and may 
be present in the inner parts of the jets: however, such winds 
cannot sustain the observed mass lo ss rates and cannot th ere- 
fore be the main engine of the jets dFerreira et alJl20()6ah . To 
date, two ac cretion powered wind models exist: i) the X-wind 
model (e.g. Shu et al. 2000h and the ii) extended disk wind 
model (e.g. iPudritz & NormanI [198611 Wardle & Koeni di ll 9931: 
iFerreira & Pelletieij|l993h . Both models are based on the same 
mechanism, the so-called magneto-rotational launching (orig- 
inally developed by [Blandford & Payne 1982 for the case of 
AGN) and only differ by the origin and configuration of the mag- 
netic field threading the disk and the size of the launching region. 
It will require higher angular resolution observations to have 
a definite ans wer regarding the process(es) at play. However, 
iFerreira et alj (2006a) have gathered indirect evidence (from jet 
rotation velocities) that appears to favor extended disk wind the- 
ories. 

If the jets are indeed accretion powered — which seems to be 
the case from our present knowledge — then the jets must affect 
the structure of the underlying region of the disk that is powering 
them. Hence, the SAD model cannot be used in this region where 
most of the angular momentum is transported away vertically, in 
the jets. In this work, we calculate the radial structure of a Jet 
Emitting Disk (JED) and compare it to that of a SAD. The paper 
is organized as follows: 
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- In Sec. |2] we briefly present the framework of the 
Magnetized Accretion Ejection Structures (MAES). This ex- 
tended disk wind model is a self-consistent description of the 
accretion disk and the jet it powers. This approach allows us 
to quantify the eff'ect of the jet on the disk. 

- Using the key parameters of the previous model, the main 
equations used to calculate the JED/SAD structure are pro- 
vided in Sec. [3] 

- The results, and in particular a comparison between the two 
(jet emitting and standard) types of accretion disks, are pre- 
sented in Sec. |4] 

- Before concluding, Sec.|5]raises a few issues and focusses on 
the possible implications of the existence of JED with regard 
to planet formation and migration. 



2. Jet Emitting Disks: The IVIAES framework 

The Magnetised Accre tion Ejection Structures model 
dPerreira & Pelletieil Il993l) has been developed so as to treat 
consistently both the accretion disk and the jet it generates. The 
idea is the same as in earlier studies of magneto-centrifugally 
launched disk winds (Blandford & Payne 1982). However, in 
the MAES, the solution starts from the midplane of the resistive 
MHD disk and evolves outwards in the ideal MHD wind/jet. 
This differs drastically from other studies where the disk was 
only treated as a boundary condition, hence forbidding any 
precise quantification of the eflect of the MHD wind on the disk. 



2.1. Properties relevant to this work 

It would be lengthy but also irrelevant to present the MAES 
model in great detail in this paper. Many papers have dealt 
with the subject, from both analytical and numerical point of 
views, and we refer the r e ader to the se papers for fu rther de- 
tails ( e.g. Ferreira 19971: Casse & Ferreira 2000ailbl: iFerreiral 
20021: iFerreira & Cassd l2004 ICasse & KeppensI l2002i l2004t 
Zanni et al. I l2007h . Instead, we give hereafter the few key ele- 
ments of the model that are important to our work. 

In a resistive MAES disk, both angular momentum and mag- 
netic field are transported using an alpha prescription. To that 
end, a local turbulent resistivity (v„, magnetic diffusivity) is sup- 
plemented to the usual turbulen t visco sity v,, (used in the stan- 
dard theorv. Sha kura & Sunvaevil973h : 



V,, = QyO-Kh^ , 

Vm = a,„VAh , 



(1) 



with Q.K the Keplerian rotation frequency, va the Alfven velocity 
and h the disk half thickness. 

Accretion solely depends on the removal of the disk angu- 
lar momentum. A relevant quantity is then the ratio A of the jet 
torque to the viscous torque. This ratio characterizes the domi- 
nant agent responsible for the extraction of angular momentum. 
For a SAD, no jet is present: angular momentum is only ra- 
dially transported away by the turbulent viscosity and A = 0. 
However, in a JED, angular momentum is vertically transported 
in the wind, along the magnetic field lines, and A > 0. It has been 
shown in previous work that steady ejection requir es A ~ 1/e, 
with e - h/r the disk aspect ratio jFerreir al l2002h . In general 
MAES solutions present large values of A: this emphasizes the 
fact that, from the moment a jet is launched, most of the en- 
ergy and angular momentum will be evacuated in the jets (see 
AppendixlAll. 



Concerning the transport of matter, the accretion rate in 
MAES solutions follows Mg oc r^, where ^ is the ejection ef- 
ficiency. One of the results of the MAES is that the latter is 
found to lie in the range [10^^, 5 x 10"'] in order to steadily pro- 
vide super- Alfvenic jets. Mass conservation in a MAES ranging 
from the inner radius r,„ to the outer radius rj reads M{rj) = 
2Mjet + M{rin) from which one gets. 



2M, 



M(rj) \ rj I Kin 



(2) 



The amount of matter being ejected thus depends on the ejection 
index, but also on the extension of the jet emitting disk. Hence, 
with the typical values r,„ - 0.04 AU, ry = 0.5 AU and ^ = 0.05, 
~ 10% of the mass will escape in the jets, in agree ment with cur- 
rent ob servational estimates (for more details see IFerreira et al.l 
l2006ah . 

The more angular momentum is removed the larger the ac- 
cretion velocity uq. The sonic Mach number in the disk = 
uq I Q^h can be rewritten under the form 



Up 



— a,,e + Iq/i — Q',,e(l + A) , 



(3) 



where a^.e denotes the effect of standard transport and 2qfi the 
specific contribution of the magnetic torque. The magnetizatiorQ 
ju - B^ljipP measures the strength of the magnetic field in the 
disk and q - figJi-h/B- is the normalized radial electric cur- 
rent density flowing at the disk midplane. This last parameter 
measures the magnetic shear as it provides an estimate of the 
toroidal magnetic field component at the disk surface, namely 

In a SAD, - a^.e, with typically a,, - 10"^ and e = h/r < 
1, so that the accretion velocity is largely subsonic. However, the 
situation in a JED is very different as MAES solutions present 
high accretion velocities with nis ~ 1. This is because steady- 
state MAES solutions are generally found c lose to equipartition 
withu e [0.1-11 and q of the order of unity ( IFerreira & PelletieJ 
11995 : lFerreira|[T997l) . In turn, q ~ I is only possible for a large 
level of turbulence, namely a,„ ~ 1: for smaller values of the 
magnetic diffusivity, the toroidal field is much larger and the 
vertical equilibrium is no more possible. As a consequence the 
torques ratio writes A = Iqiu/oye ^ where we assumed 
for simplicity ff,, - a,„. Note that if the viscous parameter a,, is 
much smaller than unity, A becomes even larger 

2.2. JED— SAD transition 

A schematic representation of the accretion structure we con- 
sider here is represented in Fig. [T] The outer parts of the disk 
have the characteristics of a Standard Accretion Disk whereas 
the inner part is occupied by a Jet Emitting Disk. The latter has 
the properties of the MAES as detailed above. The magnetiza- 
tion ji provides the criterium for the transition between the two 
types of dislQ- Let us examine this point. 

As said before, a SAD assumes the presence of a turbulent 
angular momentum transport such that the effective Reynolds 
number = ruo/v,, is of order unity. The usual assumption is 
then a comparable magnetic Reynolds number - ruo/vm, 
namely a magnetic Prandtl number of order unity (see e.g. 
Heyvaerts et al. 199^. If one considers the presence of a large 



' It is directly linked to the usual plasma beta parameter by = 2//3. 
^ Such a picture h as been also successfully applied to X-ray binaries 
IFerreira et al J2006bl . 
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no wind 








Fig. 1. Sketch of the accretion configuration suggested in this 
work. The accretion disk is constituted by a Standard Accretion 
Disk in the outer part whereas the inner part is occupied by a Jet 
Emitting Disk. See text for details. For completeness, accretion 
columns onto the central object and a possible stellar wind are 
also represented, although they are not considered in the present 
study. 



scale magnetic field threading the disc, such a small value of %„ 
translates into straight field lines, i.e. » B^. Now, the radial 
distribution B,{r) is provided by the induction equation which 
describes the interplay between advection and diffusion. Since 
there is no significant bending in a SAD, this equation writes in 
steady state 



Vm-r^ - I'lB- 
or 

with the obvious exact solution 



B^ 



(4) 



(5) 



Any large scale magnetic field will thus be naturally increasing 
towards the center in a SAD, as long as MHD prevails of course. 
Apart from triggering the MRI, the dynamical importance of this 
field is in fact measured by the disk magnetization fi{r). We thus 
need to evaluate as well the radial distribution of the total pres- 
sure P - pQ.Kh^. It can be written as 



P = 



MA 
Anrm^ 



Not all accretion disks might achieve such a transition as it 
depends on the magnetization at the disk outer edge. This outer 
boundary condition is itself reminiscent of the history of the pro- 
tostellar system, that is of both the magnetization of the parent 
cloud and the subsequent collapse. It is interesting that, indeed, 
not all young stars have detectable jets. In our picture, that would 
be explained by the lack of a JED in the in ner accretion disk. 
Using a sample of CTTS. lMenard & Duchgn e (2004) found that 
CTTS are oriented randomly with respect to the local interstel- 
lar field. This may indicate that interstellar magnetic fields play 
no strong role in enforcing the direction of the final (i.e. stellar) 
angular momentum. However, sources with strong outflows do 
have disks mostly pe rpendicular to the fi eld (i.e. jets are aligned 
to it as first found i7v lStrom et al.l ( Il986l) ). whereas sources with 
no jet detected are parallel. That could be a hint that only objects 
with disks perpendicular to the interstellar magnetic field give 
birth to magnetized central regions, namely JEDs. 

For the purpose of this paper, we will treat o as a free pa- 
rameter and compute the radial structure of JEDs. It is notewor- 
thy that rj might be observ ationally determined by measuring 
the an gular velocity in jets (iBacciotti et al.ll2002t iPesenti et al.l 
2004. 



Taking kinem atic constraints from several T Tauri jets, 
Ferreira et al.l (l2006al) found that rj would typically range be- 



(6) 



tween 0.2 to a few astronomical units. However, these values 
must be considered with caution as they might represent upper 
limits only. Furthermore, we restrict ourselves to the optically 
thick case in this preliminary study. This is justified as most 
CTTS show optically thick inner disks. Neverthel ess, one has 
to bear in mind the existence of transitional di sks (ICalvet et al.l 
120021: ID' Alessio et al.l2005l;lRatzka etaLll2007l) : the spectral en- 
ergy distributions of these disks suggest the presence of large 
optically thin inner holes, that are generally explained by invok- 
ing planet growth in these inner regions. An optically thin JED 
could be an alternative/complementary interpretation to the large 
inner gaps they display and this hypothesis deserves further in- 
vestigation. 

There are also some spectroscopic indications that gas is 
settled in th e innermost hot d isk regions where dust has been 
sublimated dNaiita et al.ll2007l) . However, we will not consider 
line emission in this work as these inner regions are probably 
strongly affected by the stellar magnetosphere. We will therefore 
use for the disk inner radius r„, an arbitrary but representative 
value of the disk truncation radius, r,>, - 0.04 AU. 



where both the accretion rate Mg and nis - a^.e have been as- 
sumed constant Using Eq. (|5]l, we obtain 



HOC r ^ with 



6 = 2%n-5l2 



(7) 3. Calculation of the disk structure 



Strictly speaking, the magnetization will increase towards the 
center whenever %„ > 5/4. Thus, unless the magnetic Prandtl 
number is significantly smaller than unity, this is a condition that 
is certainly automatically satisfied in SADs. However, making 
a transition to a JED requires a large value of /y, namely fi ~ 
1; a JED will therefore exist only below a transition radius rj 
such that i-iirj) ^ 1. Below this radius, the radial dependency of 
the vertical magnetic field is different than that given by Eq.©. 
Indeed, a JED requires a field close to equipartition throughout 
all its extent so that one gets B^ cc y-^^'^+^l^ (^ is the ejection 
efficiency). 

' Note that e is only slowly varying with the radius in a SAD around 
a protostar so that /? oc r is a good approximation. 



We are interested in deriving the radial structure of an accretion 
disk, be it a SAD or JED. The advantage of the approach used in 
this work is that both types of disks are described using the same 
formalism: the only difference lies in the values of the MAES 
parameter A, thus nis- We make the assumption of a geometri- 
cally thin, optically thick steady-state disk, rotating at Keplerian 
velocity, i.e. with a frequency D.k - ^Jgm'jT^. We also assume 
that the accretion rate Mg does not depend on the distance to 
the central object. This is justified since soluti ons to the MAES 
proble m give Mg oc r^, with^ ~ lO^-'-SxlO"' dCasse & Ferreira! 
I2000bl) . We also assume the gas and the dust to be well coupled, 
at the same temperature, and that the mixture behaves as an ideal 
gas. 
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3.1. Basic equations 

We restrict ourselves to the simple case of a steady state ac- 
cretion disk. The calculation of the disk structure relies on the 
equality of the cooling and heating terms. The latter can be writ- 
ten in a generic form as 



<2^ = / X 



(8) 



where / represents the fraction of gravitational potential energy 
that is used to heat the gas. For this work, the gas is only locally 
heated by viscous effects and any other source of heating, such 
as irradiation from the staiQ or dust-gas collisional heating (oc- 
curring when Tdust ^ Tgas) are discarded. In the MAES context, 
it has been shown (see AppendixlAli that 



1 



1 + A 



(9) 



- For a SAD, A = and /sad = 1 : in that case, the viscosity 
is responsible for converting all the mechanical energy into 
heat. 

- For a steady-state JED, however, A ~ e ' » 1 and/jEo 1: 
as mentioned earlier, when a jet is present, most of the disk 
energy leaves in the jet as an MHD Poynting flux. Only a 
small fraction of that energy remains in the disk to be con- 
verted into heat and radiated away. 

As for the cooling, we assume that the disk radiates like a 
black body with an effective temperature Tg^, which leads to 
(Hubeny 199(J 



(10) 



where the optical depth r a; Kp^h links the effective temper- 
ature Teff to the mid-plane temperature Tq via the opacity of 
the gas K. This link between effective and central temperatures 
hold s only if energy t ransport is done vertically by photon diffu- 
sion. [O^Aleiiioieral] ([T998i|) showed that radiative engery trans- 
port was indeed the dominant mechanism in SAD, with respect 
to convection or turbulent transport. However, the level of tur- 
bulence required in a JED {um ~ 1) is much larger than that 
usually assumed in a SAD (a^ ~ 10"^) so that the results of 
iD'Alessio et al.l (Il998l) may not apply here. Nevertheless, with- 
out any reUable expression for a turbulent energy flux, we restrict 
ourselves to the usual approximation. For the opacity, we adopt 
the standard k = kPqTq form, where k, a and b have to be ad- 
justed regarding th e dominant coolant present in the gas. We use 
the values given in iBell&LinI (11994 that includes eight opac- 
ity regimes, from dust dominated cooling to electron scattering 
cooling (although this regime is not reached by the tempera- 
tures at play in YSO disks). 

The Bell & Lin opacity adjustments can be used to derive the 
disk structure analytically, which is not the case with tabulated 
results. The value of the opacity of the gas at a given temper- 
ature and density is a fundamental quantity regarding the cal- 
culation of an accretion disk structure and should be evaluated 
carefully. Many authors have numericall y calculated Rosseland 
mean opacities, us ing different models (Alexander & Ferguson' 
Il994t iHenning & Stognienko 1996; Semenov et al. .2003.) and 



* This hypothesis is further discussed in Sec. 15. II 

^ In Appendix|B] we give the general expressions of the radial struc- 



the Bell & Lin prescription deviates from more refined ap- 
proaches ar ound 1500-1800 K, where it significantly under- 
estimates K d Semenov et al. I I2OOI . Thus, to be consistent, we 
stopped our disk calculations whenever the central temperature 
reached this value. For all JED models, the inner disk edge ri,, 
is achieved at lower temperatures (see Fig. |2|. However, SAD 
models are affected and this is the reason why, in that case, the 
disk inner radius is actually l arger than the chosen r,„ (see also 
iPapaloizou & Terquemlll999h . Hence, more refined SAD calcu- 
lations should be performed at these inner radii, probably taking 
into account dust sublimation and line emission from the gas 
component. This is however beyond the scope of the present pa- 
per 

3.2. Radial structure of fhe accretion disk 

A standard w ay of writing the vertical hydro static equilibrium of 
a disk is (e.g.. IPapaloizou & Terquemll 1 9991) 



dP 



Integrated over the thin disk, the hydrostatic equilibrium defines 
the scale height of the disk h, via the aspect ratio e - h/r: 

2 ^0 



where quantities with the subscript 0, refer to mid-plane values. 
Using the perfect gas equation of state, the previous equation 
reads 



kBTo 



(11) 



where ks is the Boltzmann constant, fi the mean molecular 
weight of the gas and nip the proton mass. 

Defining t/o to be the inward radial velocity of the flow, one 
can simply calculate the density as 



1 



po = 



(12) 



AnrhuQ AnQ.K^' nise^ 

where m, is the sonic Mach number, defined by Eq. (O in the 
MAES model. 

Equaling the heating and cooling terms 



(13) 



and using Eq. ( fTTT i and Eq. ( fTSl i it is possible to express all the 
thermodynamical quantities (in particular the midplane temper- 
ature To and the disk surface density Z = 2p/!) as a function of 
the radius and of the parameters of the problem, namely under 
the generic form 



X(r) oc J?(^.(M*, Mfl, nis, k, a, b) r"' 



(14) 



ture of a JED, keeping k, a and b as free parameters. 



where and 6x relates to the quantity X and depends on the 
problem parameters. However, the explicit forms of these quan- 
tities can be very lengthy and are postponed to Appendix IB] for 
the sake of legibility. Note that these expressions could be di- 
rectly used by those wishing to quickly derive the structure of a 
JED for a given set of stellar and opacity parameters. 

Concerning the magnetic field and as stressed in ^2.11 the 
necessary condition for launching a self-collimated jet from a 
Keplerian accretion disk is the presence of a large scale vertical 
magnetic field close to equipartition (// ~ 1 . iFerreira & Pelletierl 



Combet & Ferreira: Jets and accretion disk structure 



5 



Il995h . This condition allows us to estimate the strength of the 
magnetic field 



10^ 



IMA^'U Ma r Y^i^^^n f::::;::::::::;:::::::::::::::::::::::::;::::r' h ^>ii 

We remind that for the specific case of our calculation, the eiec- 1 -++++++++++++++++++++++++++++++ 

A ^ _Y --++++ +++++++++ ++++++++++++++ 

tion efliciency ^ = 0. This equation along with the set provided & 10 ?+++♦++++++++♦*+++++++*+ 

in AppendixlBl(rn. e and 2) completely determine the disk struc- • S 

ture. :-;;n";;;;;r" 
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4. Results 



The main features of jet emitting disk can now be derived from 10 



10' lo" lo' 10" 



the expression collected in Appendix |B] In Sec. 14.11 the ra- r (AU) 

dial variations of the key quantities of the JED are presented. 

Comparison to the standard case is also drawn. Sec. O illus- Fig. 3. Region of optically thick JED in the (r, M) plane. Crosses 
trates the eff^ect of a jet emitting region on the disk on the spectral represent the validity range of our calculation, 
energy distribution. 



4.1. JED: main characteristics 

In Fig. |2] are plotted in solid lines the radial dependences of the 
central temperature (upper-left), aspect ratio (lower-left), surface 
density (upper-right) and magnetic field intensity (lower-right) 
of the JED, as given by Eq. dBJl i, (EH), dES and The 
quantities — ^but for the magnetic field — appear as broken power- 
laws, each segment corresponding to a given opacity regime 
(k, a, b). The transition radius between two regimes of the Bell 
& Lin opacity, say ; and / + 1, are obtained by Ki{r) — Ki+i{r) and 
solving for r. 

The disk quantities have been plotted for different mass ac- 
cretion rates, from 10"^ to 10"^ yr"'. In agreement with 
intuition, the higher the accretion rate, the higher the tempera- 
ture, surface density and disk aspect ratio. Note also that more 
opacity regimes need to be taken into account for the highest 
accretion rates as a result of the increasing the temperature. We 
have truncated each curve at the outer radius for which the gas 
becomes locally optically thin to its own radiation and where 
our description of the radiation in Eq. ( fTOb stops being valid. 
But since accretion disks of CTTS are optically thick, this ra- 
dius marks an upper limit for rj. Figure [3] is another represen- 
tation of the latter point and gives in the (r, M„) plane the re- 
gion where the optical depth t is greater than unity in a JED. 
There are not many sources where both the accretion rate and 
rj have been estimated. This is however the case for DG Tau 
{Mg ~ 2 X 10''' M e yx-\ iBacciotti et all |2002|) for which 
iFerreira et al.l ( l2006al) estimated rj ~ 1.5-4.5 AU: these val- 
ues fall into the validity range of our calculation given in Fig.[3j 

Coming back to Fig. |2l all quantities (but for the magnetic 
field) are compared to the SAD case (dashed curve). For the 
latter, we have plotted only the highest accretion rate (and use 
a,, = 10'^) to prevent overlapping curves. Also not shown here, 
our SA D models are in agreement w ith the disk structure de- 
rived in lPapaloizou & TerquemI ( Il999l) . For the SAD, curves for 
smaller values of the accretion rate and other a,, are available in 
that paper. 

The radial structure of a JED is found very different than 
that of its equivalent SAD: at a given accretion rate, the JED 
is cooler (Fig. |2^), lighter (Fig. [2}5) and thinner (Fig. [2j;) than 
the SAD. In particular, at a given radius, the surface density can 
vary by ~ two orders of magnitude. This implies a large density 
jump at the transition radius o, that remains whatever the chosen 



value of the turbulent parameter a,.. Implications of this point are 
further discussed in ^5.21 

Turning to the magnetic field, one might object that fields of 
the strength given by Eq. ( fTSl l and illustrated in Fig. [2}l are im- 
possible in accreting systems. However, the value of this mag- 
netic field is actually far smaller than the one estimated from the 
interstellar ma gnetic field assum i ng either ideal MHD oc n 



or B cc ni/^ ( iHeiles et al.lll993l: iBasu & Mouschoviaslll994l) . 



Indeed, if one takes the fiducial values n ~ 1 cm'^ and B ~ 4 
observed within dense clouds and use the law B ccn^l'^ jCrutchej 
1 19991) . we get a magnetic field at 1 AU ranging from 10 to lO"' G 
(depending on the density)! Thus, the main problem is to get 
rid off" the magnetic field during the infalling stage. This is- 
sue is still under debate. It seems nevertheless straightforward 
to build up accretion disks threaded by a large scale magnetic 
field of a large amplitude from 3D collapse calculations (see e.g. 
iBaner jee & Pudritz 2006). 

The observation of magnetic fields is a very difficult task and 
no measurements have been obtained so fa r in a disk known 
to drive a jet. However, iDonati et al.l (l2005l) have managed to 
measure the magnetic field strength in the disk of FUOr, a very 
strong accreting object with Mg ~ lO'^M© yr"'. Using the 
spectro-polarimeter ESPadOnS, they found a value of B ~ 1 kG 
at 0.05 AU from the star At this distance and for such an accre- 
tion rate, Eq. ( flSl l gives B ~ 100 G. The magnetic field present 
at the inner disk of FUOr is actually higher than equipartition: 
according to the MAES theory, no steady state self-collimated 
jet can be launched from this disk, which is indeed consistent 
with observations. If one takes this observation at face value, 
then one possible explanation is that the disk magnetic flux has 
been advected and compressed towards the star by the sudden 
rise of disk accretion rate. This conjecture should deserve fur- 
ther investigation. Note also that the presence of a strong (i.e. 
larger than equipartition) vertical magnetic field in an accre- 
tion disc triggers no n-axisymmetric instabilities and spiral waves 
(iTagger et al.ll992l) . If these waves bounce back at the inner disk 
boundary then a standing spiral pattern is formed and leads to the 
formation of a magnetized vortex (a MHD Rossby wave) local- 
ized at the radius where the keplerian rotation coincides with the 
wave frequency ([Tag ger & Pellat 1999; Caunt & Tagger 2001]; 



Stehle & Spruilj|2q0ll) . Although the field structure derived by 
Donati et al. shows a remarkably high degree of axisym- 
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Fig. 2. Radial variation of key disk quantities. Upper-left: mid-plane temperature. Upper-right: surface density. Lower-left: aspect 
ratio of the disk. Lower-right: large scale vertical magnetic field intensity. Line thickness represents the increasing accretion rates 
considered for the JED: M„ - 10"^ , 10"^ , 10"'' and 10"^ Mo yr"'. However, the SAD case is computed for Ma - 10"^ Mq yr"' 
and ay - 0.01 only in order to prevent overlapping curves. 



metry, one might indeed expect non axisymmetric perturbations 
in strongly magnetized disks. However, for the physical condi- 
tions envisioned in JEDs, namely a field close to but smaller 
than equipartition, such instabilities are rather weak (see how- 
ever lKeppens et al.lllool for MHD instabilities at ju ~ 1). 

4.2. Spectral Energy Distributions 

The physical properties of a jet emitting disk are rather differ- 
ent from that of a standard disk. In consequence, their radia- 
tive properties should also differ and are investigated hereafter. 
Spectral energy distributions (SED) are one of the main disk di- 
agnosis for comparison with observations. If the radial structure 
of the disk does not depend too drastically on the irradiation of 
the central object (see ^5.1 1, its radiative properties certainly do 
(see iDullemond et al.ll2007 for a review). In particular, the sim- 
plest approach considers that the stellar illumination can create 
a super-heated laye r at the surface of the disk which, in turn, af- 
fects the disk SED (IChiang & Goldreichlll997l:lDullemond et al.l 
I2OOI . The success of this model comes from both its simplic- 
ity (compared to that of a full radiative transfer treatment) and 
ability to reproduce observations of TTauris and Herbig Ae/Be 
stars. 

However our purpose here is to compute the SED to illustrate 
a possible effect of the JED rather than reproduce specific obser- 
vations. For that reason, the irradiation of the central star on the 
disk will be included the very crude fashion detailed below. . 



We will consider two components to the irradiating flux: the 
stellar radiation and the UV accretion luminosity. Indeed, in the 
now widely accepted magnetospheric accretion model, matter is 
channelled a long the stellar magn etic dipole from the disk onto 
the star (see iBouvier et alj|2007l and references therein) where 
it creates the observed hot spots/rings. The UV radiation from 
these accretion shocks can then illuminate and heat the disk. In 
the following, Luv will refer to the UV luminosi ty of one on the 
two ho t rings, and to the stellar component. iGuUbring et al.l 
(Il998h estimated from a sample of TTauri stars that 



2L, 



3.5 ' 



where the total accretion luminosity scales with the accretion 
mass rate as 



Lncc ~ „ I ^ 



The stellar luminosity is simply given by 



(16) 



(17) 



where and are respectively the star radius and tempera- 
ture. The irradiation flux Q^^ is then obtained by 



Qn 



Li, +L 



UV 



Anr^ 



cos ifj , 



(18) 
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with ip the angle between the incident radiation and the normal 
to the disk surface. For the sake of simplicity, we will assume a 
fixed incident angle, cos ijj - 0.05 (i.e., incident angle ~ 3 deg.): 
this is of course not consistent with our disk model, as /i/r 9^ est 
(Fig. [2];), but has the benefit of giving a quick estimate of the 
irradiation component. 

With these ingredients, the disk effective temperature re- 
quired by SED calculation is obtained by 

crT'^f, = e+ + gin- . (19) 

The viscous heating flux is given by Eq. ([8]l — we remind here 
that this flux is much smaller in a JED than in a SAD. 

To actually compute the SED, we consider the disk configu- 
ration presented in Fig.[Tl where a JED is present from rin to rj 
and a S AD from rj to rpm. W e choose arbitrarily ry = 0.5 AU 
(see e.g. iFerreira et al.ll2006ah . For illustrative purpose. Fig. 
gives the disk height scale as a function of the radius. The thick 
solid line corresponds to the disk studied here, with a transition 
in rj. The SED obtained from this structure is then compared 
to the case where a standard disk is present on the entire spatial 
domain [^n, rom]- Using the temperature derived in Eq. (T% . the 
disk SED is computed for a one solar mass star and using the 
simplest geometry, where the disk is viewed pole-on, 

rBv(reff) dr 



= Sn'v{^j^' rB,{TiED)dr+ jj rB,(TsAD) drj . (20) 

In this expression, Planck's law By depends of the radius through 
the disk effective temperature T. The contribution of the central 
object to the SED is also considered, using - 4000 K and 
Ri, = 2.5 Rq, and assuming it radiates like a spherical black- 
body; 

vLT' = 47T^RlvBy(T^) . (21) 

The spectral energy distribution is plotted in Fig.|33, from in- 
frared to millimeter wavelengths. The dotted and dashed curves 
correspond respectively to the stellar and disk contributions. 
The total SED is plotted in thick solid curves, for each disk 
configuration — SAD alone or JED+SAD. The effect of a jet 
emitting region in the disk is characterized by a smaller flux and 
a redshifted disk SED. Note that Kuncic & Bicknell (20071) very 
recently reached a similar conclusion more specific to the case 
of AGN. However, although they used similar ingredients as in 
the present study, their outflow model was purely phenomeno- 
logical and did not rely on a physically consistent underlying 
model (such as the MAES model) for the quantification of the 
jet torque. This effect increases with the size of the transition 
radius rj: this is shown in Fig. |4]:, where the discrepancy be- 
tween the standard SED and the JED+SAD becomes larger for 
larger rj. This is in agreement with intuition as it was shown in 
the previous section that JED were lighter and cooler than their 
equivalent SAD. 

The last effect on the SED put into light in this section is 
shown in Fig. |4}l. For a fixed transition radius rj - 0.3 AU, 
the SAD and JED+SAD spectral energy distributions are com- 
puted for several mass accretion rates, namely 10"^, 10"* and 
lO^'^Mo yr-'. For a given line pattern, the thickest line corre- 
sponds to the standard disk alone. For the highest accretion rate, 
the standard disk SED completely overwhelms the stellar com- 
ponent which explains why the thick dashed curve does not join 
the others at small wavelength. 



Two interesting effects are put into light in this figure. First 
of all it is found that the discrepancy between the standard and 
JED+SAD energy distributions is higher for higher accretion 
rates: at IQ^^Mq yr"' (and for that particular rj), the standard 
SED is at most ten times larger than the JED+SAD one, whereas 
this value falls to 2.5 at IO^^Mq yr"' . The second effect, linked to 
the first, concerns the location of this maximum discrepency. It is 
found around l.S/im at 10"*' Mq yr"' and ~ 3/im at lO^^M© yr"'. 
These two points may be relevant regarding observational efforts 
on the matter However, let us stress one more time that these 
SED are presented to enlight a potential effect only. Further work 
is required to reach a real quantification of this effect when tak- 
ing properly into account the central object heating of the disk 
surface. 



5. Discussion 

5.1. Illumination from the central star 

In the previous section, we included irradiation of the central ob- 
ject in the calculation of the SED as it is known to have a strong 
influence on the latter However its influence onto the disk struc- 
ture is less well established. For that reason, but also for sake of 
simpUcity, the irradiation from the central star was not taken into 
account in our calculation of the disk structure. Nevertheless, 
as it decreases in r^^ compared to the r^^ of the viscous heat- 
ing Eq. (ISj, it is expected to become dominant beyond a certain 
distance from the c entral object, geometrical effects put aside. 
iDubus et~an (1 19991) emphasized that the appropriate criterium 
for the disk structure to be dominated by viscous heating reads 



(22) 



where Qm is the irradiation flux coming from the star (see 
also iRafikov & De Collell2006l) . The total optical depth of the 
disk Ttot appears as the irradiation mainly affects the surface of 
the disk whereas the viscous heating is assumed to be present on 
the entire disk thickness. 

The irradiation flux has been determined in the previous sec- 
tion, by Eq. ( fTSl ) and can readily be used with Eq. ([8]l to check 
the validity of our disk structure calculation given by the above 
criterium. However, for this calculation, we will make no as- 
sumption on the incident angle for the radiation but compute it 
consistently given the disk structure we previously derived. The 
source of the radiation is considered point-like. This is a reason- 
able assumption for the UV hot spot located at the end of the 
accretion column. This is however questionable for the stellar 
component: the star is an extended source, at least for the most 
inner radii. Again, for simplicity, we will consider the two com- 
ponents (stellar and UV) as coming from one unique location. 
By changing the altitude of the spot, we should be quite conser- 
vative in our conclusions. 

From geometrical considerations, one gets 



cos if/ 



d\nh ^ ^ Zsh 
dlnr r 



(23) 



In the previous expression, all quantities have their usual mean- 
ing and Zsh i s the light spot alt itude. If Zsh = 0, the standard 
expression of iFrank et alj (|2002|) is recovered. 

Within the JED, we found that when z^h equals zero the 
illumination could safely be discarded for any accretion rate 
Ma ^ 3 X IO^^Mq yr"' (at any radius). The situation changes 
when the hot spot is located at a higher altitude as it can strike 
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Fig. 4. Upper-left: disk scale height as a function of the distance to the central object. The JED/SAD transition radius is ry = 0.5 AU. 
Upper-right: pole-on viewed spectral energy of the JED+SAD disk of (a). Dashed lines correspond to the SAD (thin red) and 
JED+SAD (thick black) configurations. A 4000 K stellar blackbody is included (dotted line). The sum of the stellar and disk SED 
is plotted in thick solid curves (red and black). Lower-left: variation of the total JED+SAD energy distribution with increasing 
transition radius. Lower-right: influence of the accretion rate on the SED for a fixed transition radius. Thick lines correspond to the 
SAD alone and thin ones to the JED+SAD case. 



the disk more directly in the inner regions. For Zsh - 0.77?*, i.e. 
an accretion column starting at the star at a colatitude of 45°, we 
found that Eq. (I22]i was not satisfied for M,, < 7 x IO^^Mq yr 
To obtain these thresholds, we required Eq. (l22l i to be true on 
the entire radius range. Note that these accretion rates could be 
lowered by some extent if one only requires the condition to be 
realised on most of the radius range: it could be argued that a 
narrow irradiation dominated region will not change much in 
the overall structure of the disk. Basically, illumination of the 
JED can be neglected at any radius for M > a few lO^^M© yr"'. 
This threshold is consistent with accretion rates of sources driv- 
ing powerful (detected) jets. Nevertheless, it would be of interest 
to extend this work to the lower accretion rates by including ir- 
radiation in the heating terrrQ. From these considerations, we are 
left confident with the validity of our very simplified analytical 
description, at least for M > a few IO^^Mq yr"'. 

5.2. Implications for planet formation and migration 

The initial conditions for planet formation depend on the disk 
physical properties. It has been shown in the previous section 
that the latter are very different from a JED to a SAD. We dis- 



* This can only be done iteratively (hence, not analytically) as the 
knowledge of the disk structure is required for the determination of 
cosi^. 



cuss, hereafter, two of the possible consequences that an accre- 
tion structure such as that of Fig. [T] may have on planet forma- 
tion. 

iGammid (Il996h was the first to present the idea of layered 
accretion disks where the upper part of the disk is ionized — via 
collisions, cosmic rays or X-rays — and some embedded inner 
part stays neutral. The latter, termed dead zone is then decou- 
pled to the magnetic field, hence to the MRI induced turbulence, 
and remains quiescent. Angular momentum is very poorly trans- 
ported outward in the dead zone and no accretion occurs in this 
region of the disk. 

We have shown in the previous section that jet emitting disks 
were both thinner and lighter than standard disks. As a conse- 
quence, their are likely to be more ionized than SAD as X-ray 
radiation from the central star and cosmic rays should deeperly 
penetrate them. The calculation of the actual ionization degree 
of JED is postponed to a forthcoming study, which should give 
a definite answer to the possibility of a dead zone in a JED. 
However, it is of interest to briefly mention a few issues related 
to that matter. 

Both low and high mass planet formation start with the 
growth of planetesimals and the necessity of dust settling and 
agglomeration. The effect of turbulence on agglomeration is still 
being debated. Indeed, it might help to for m planetesimals by 
trapp ing dust within turbulent vortices (e.g. lBarge & Sommerial 
Il995l) . On the other hand, recent numerical simulations of a 
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fully turbulent disk showed that only the larger grains (with 
typical sizes from 1 to 10 cm) do settle towards the mid- 
plane whereas the smaller on es stay in suspension in the disk 
dFromang & Papaloizoull2006l) . The thickness of the dust sub- 
disk was also found to be rather large (~ 0.2/1). These authors 
also showed that the dust subdisk is thinner (more compact) 
when a dead zone is taken into account in the simulations. Along 
the same line, iCie sla ( 2007) stressed that the efficiency of dust 
settling and coagulation in the presence of a dead zone was en- 
hanced compared to that of a ionized (i.e., turbulent) disk. From 
those considerations, layered accretion disks appear to be a priv- 
ileged location for grain growth and planetesimal formation. If 
JEDs are indeed fully ionized hence turbulent, then it is unlikely 
that they can host the earliest stages of planet formation. 



On the other hand, the maximum value of the JED-SAD tran- 
sition radius o, as determined by the transition to the optically 
thin JED regime, is rather small (see Fig. [3). It reaches ~ 10 AU 
only for high accretion rates, found only in embedded sources 
(Class or 1). For accretion rates of some 10"^ yr"', more 
typical of CTTS, rj cannot exceed 0.1-0.3 AU. As the JED has 
a rather low column density 2 (see Fig. [2j'), it provides very 
unideal conditions for the core accretion process leading to mas- 
sive planet formation. 



Planet migration in a gaseo us disk surrounding a central star 
(iGoldreich & Tremaind [19801) has been a long standing issue 
regarding planetary formation theories. The transfer of angu- 
lar momentum between the planet and the gaseous disk gen- 
erally results in the inward motion of the planet. This general 
mechanism affects both low (ty pe 1 migration) and high mass 
planets (type 11 migration) — see (lArmitagel2007l) for a review — 
and is generally invoked to explain the existence of hot Jupiters. 
However, the existence of giant planets far from their central 
star (a few AU) is a real issue. Indeed, type 1 migration is known 
to be a very fast process and should not allow enough time for 
core accretion (build up of a massive solid core), i.e. for the ex- 
istence of giant planets at such distances. Hence, finding a way 
to slow down or halt type 1 migration is of importance if one is 
to understand planetary system formation. 



In that respect, iMasset et al.l (l2006l) studied the effect of a 
sudden radial surface density decrease (going inwards) in the 
disk on type 1 planet migration. These authors found that such a 
density jump could indeed trap a protoplanet at the location of 
the transition by a balance between the corotation and Lindbald 
torques. As shown previously, the surface density of a SAD is 
always larger than that of a JED for a given accretion rate (see 
Fig. lit). Hence, a transition from an outer SAD to an inner JED, 
as described in Sec. 12. 2! would naturally provide the surface den- 
sity jump necessary to a planet trap. It is therefore likely that 
planetesimals should first form in the outer, denser SAD and 
then migrate towards the center to be halted at the JED-SAD 
transition. 



One must however be cautious as in lMasset et al.l (l2006h the 
physics of the JED and, in particular the presence of the large 
scale magnetic field, was not taken into account. It would be 
useful, as a second step, to include this critical ingredient and 
determine precisely how will the JED be itself affected by the 
presence of a protoplanet. 



5.3. The Ma-Mi, relationship in the JED model 

One of the striking observational properties of T Tauri stars and 
brown dwarfs is the apparent steep correlation between the mass 
of the central object and the disk accretion rate, namely 



M^ocMZ 



(24) 



with a lying between 1 and 2 (iMuzerpUe et aP l2003 | : 



' Calyet et al.l2004l:lAiexander & Armita ge 2006': lHartmann et al.l 
120061) . Though seemingly very robust observationally - four 
orders of magnitude in accretion rates and two i n mass-, the 
physical origin of this relation is still debated. iPadoan et al.l 
(i2005i) evoked the possibility of a star-disk accreting from a 
large scale envelope at the Bondi-Hoyle rate, giving precisely 
Ma oc M^ . However, this scenario does not take into account the 
angular momentum of the infalling gas. In an other approach, 
iDullemond et al.l (l2006l) showed how taking into account the im- 
print of the physical properties of the parent core onto the star- 
disk system lead to Ma k M|-^. lHartmann et al.l (l2006b explored 
several mechanisms but concluded on a word of caution: these 
authors noted that, if the bulk of the data is well represented by 
M„ oc Ml, the TTauri with th e highest accretion rate s appear to 
have Ma <^ M^,. In that sense. iHartmann et al] (|2006|) stated that 
it may not be relevant to look for a universal Mg-Mi, relation- 
ship, as different mass regimes may have different dominating 
processes. Nevertheless, we believe it is of interest to see what 
Mg-Mi, relationship is predicted by the JED model. This is done 
below, using very simple arguments. 

The disk accretion rate M„, as derived from observations, is 
actually indicative of the mass flow at the inner disk regions, 
hence representative of the JED rather than of the outer SAD. 
But in a JED accretion is achieved through angular momentum 
removal in jets. This is made possible only because of the pres- 
ence of a large scale disk magnetic field Bdisk- As argued in 
Sect. 2.2 this field is probably the parent cloud core magnetic 
field Z?core that has been advected and concentrated by the in- 
falling material during the collapse. Such a collapse is induced 
when the ratio of the magnetic flux to the total mass ^/M^ore 
reaches a critical value. This translates into 



Brnre 0^ Mn, 



(25) 



as an initial condition. Now, assuming that the ratio O/M re- 
mains constant (or varies only slowly) during the collapse, one 
gets 



(26) 



since most of the infalling cloud core mass ends up in the central 
star (ejection will not change the relation oc Mcore)- Note that 
another way to obtain a relation of this kind is by writing Bdisk 
Bqok^ where the value of the exponent rj should be provided by 
full 3D collapse calculations. However, we expect it to be close 
to unity which would lead to Eq. (l26l l as well. 

The disk magnetic field required in a JED is given by 
Eq. (ITSl l and leads to 



(27) 



Inserting this relation into Eq. ( |26] l provides then the Ma-M^. 
relationship expected in JEDs, namely 



Ma oc mI'^ 



(28) 



If our crude derivation gives a result consistent with the ob- 
servations, one must however remain cautious: for a given stellar 
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mass M j, there is observat i onnaly a large dispersion in Ma- This 
has lead lClarke & Pringld ( |2006|) to point out a possible incom- 
pleteness of the samples used, that may have lead to such an 
apparent correlation. 



6. Conclusions 

There has been growing observational evidence that the bipolar 
ejections of matter occurring during the formation of a star are 
directly linked to the accretion process. From that point of view, 
jets should have an influence on the structure of the region of the 
disk that is launching them. 

Using the framework of the Magnetised Accretion Ejection 
Structures, we derived the radial structure of such a Jet Emitting 
Disk and compared it to the equivalent Standard Accretion Disk. 
We found that a JED is cooler, thinner, and lighter than a SAD 
at the same accretion rate. Invoking the idea of a radial transi- 
tion from an inner JED to an outer SAD, we have also shown 
in a very crude approach that such a transition has an eff'ect on 
the Spectral Energy Distribution: the flux coming from a JED is 
smaller and redshifted compared to the SAD. A more detailed 
study is however needed to provide realistic observational pre- 
dictions. 

Several implications on planet formation have also been 
drawn from the existence of JEDs. On the one hand, and al- 
though this needs confirmation, dead zones may not exist in 
JEDs as their small thickness and density will favor ionization. 
Recent numerical works show that dust settling could be made 
more difficult in the absence of a dead zone. This questions the 
possibility for a JED to host the earliest stages of planet forma- 
tion. On the other hand, the surface density jump occurring in a 
JED/SAD transition could serve as a planet trap and halt type 1 
migration. 

Protostellar jets have been observed and theoretically studied 
by the magneto-centrifugal approach (among others) for almost 
thirty years. Nevertheless, little attention has been paid to their 
potential effects on the disk structure. In this paper, we stress that 
these effects indeed exist and further work is to be undertaken for 
better quantification and observational perspective. 
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Appendix A: l\1AES global energy budget 

The global energy budget of a magnetized accretion-ejection 
structure writes (iFerreira & Pelletier, il993i: iCasse & Ferreiral 
I2000b^ 



— 2Pjet + 2Pr; 



(A.l) 



where the accretion power, namely the power released by the 
accreting flow in the disk 



2n„ 



rjn 
rj 



is obtained by computing the difference between the mechanical 
power that comes in at ry and goes out at r„, . Note that Maj is a 
constant and must be understood here as the accretion rate feed- 
ing the JED at ry (within the JED, the disk accretion rates varies 



as Ma °^ 1^)- All available power P^cc is thus shared between 
radiative losses at the disk surfaces Piad and jet power fjet. This 
last term marks obviously the difference with the global energy 
budget of a SAD. In the latter, all accretion power is released as 
radiation whereas JEDs also power jets. This last contribution to 
the budget is actually the easiest to compute as it is the flux of 
energy that leaves the two disk surfaces namely 



^jet = J" [Smhd + p£up] ■ dS 



(A.3) 



where £ = m^/2 + <l>c + H is the non-magnetic specific en- 
ergy (Og is the gravitational potential and H the enthalpy), 
^MHD - -^FfB^Bp/jJo is the MHD Poynting vector {Qf ^ Qk 
is the angular velocity of the magnetic field line) and dS = dSn 
with dS the elementary disk surface and n the unit vector normal 
to the disk surface. For magnetically driven jets from keplerian 
accretion disks, the dominant term is the magnetic contribution 
due to the MHD Poynting vector. Thus, the energy flux carried 
away by the jets depends directly on the toroidal field at the disk 
surface B\ and writes 



2P 



jet 



A 



(A.4) 



where A~l/e»lis the ratio of the jet to viscous torques 
and q the magnetic shear parameter. Most MAES solutions were 
found with |Z?^| ^ qB., which translates into Pace ^ 2Pjet. Thus, 
most available power feeds the jets and energy conservation tells 
that only a fraction of it, namely 



2P, 



rad 



1 



Pace 1 + A 

is dissipated within the JED. This justifies the use of Eq. (|9]l. 
Note however that some solutions were found with \B'^\lqB. as 
low as 0.5, which shows that JEDs may also, under certain cir- 
cumstances, produce some luminosity. But as a first step to illus- 
trate the effect of JEDs in the center of YSO accretion disks, we 
have disregarded these peculiar solutions (see also discussion in 
.Ferreira et al..,2006b) . 

Appendix B: JED key quantities 

In this appendix, we give the detailed analytical expressions of 
the main disk quantities as calculated in the MAES context. We 
only focus on the radial dependence of each quantity. In these ex- 
pressions, the Rosseland mean opacity is taken under the form 
K - Kp"T'', where the parameters it, a and b are to be determined 
depending on the considered opacity regime. In the following, 
cr is the Stefan-Boltzmann constant, and Ma respectively 
represent th e mass of the central object a nd the mass accretion 
rate, fl - 2 (iPap aloizou & Terquemlll999h is the mean molecu- 
lar weight of the gas, nip the proton mass and ks the Boltzmann 
constant. The parameter nis is related to the MAES model (see 
Sec. |2]l and is found to be close to unity in most solutions of 
the problem. Using Eqs. (fTTT i and (fTST i to rewrite Eq. ( fT3T l from 
Eqs. ([8]l and (fTOj i, the radial dependence of the mid-plane tem- 
perature is expressed as 



(A.2) To{r) = r 



-r[i+5(fl+i)/2] 



(B.l) 



where F = l/{4- + a - b) and is dependent of the parameters 
of the problem. 



\^ncr/ \4nms 



1 



r(fl+i) 



(GM^y 



r(fl+l)/2^r(a+2) 



kB I 
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Plugging Eq. dB.ll i into Eq. (fTTT i. the disk aspect ratio e = /z/r is 
then directly obtained under the form 



(B.2) 



-{ 



r/2 



1 



— 1 , 

SttctJ \4mn 



The disk mass density can be calculated in a similar fashion, 
(using Eqs. (IB. 21 ) and ( fT2b together), but we rather give the ex- 
pression of the surface density as it appears to be a more relevant 
quantity regarding accretion disk studies. The expression of S(r) 

X+h 
1^ p{r,z)dz ~ 2poh 

and reads 



with 



2p(i(r)e(r)r 
(B.3) 
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